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1.  Foreword 

The  Sun  is  an  abundant  source  of  energy  capable  of  meeting  all  our  energy  needs  if  properly 
harvested.  However,  efficient  capture,  storage  and  transport  of  energy  from  sunlight  still  remain  a 
challenge.  One  of  the  ways  to  utilize  solar  energy  is  conversion  of  sunlight  to  electricity  via  dye-sensitized 
solar  cells  (DSSCs)  or  to  chemical  fuels  via  photocatalytic  synthetic  cells.  These  systems  are  often 
designed  around  a  single  photoactive  molecule  (a  chromophore)  or  a  molecular  array  anchored  to  a 
semiconductor.  The  conversion  of  sunlight  to  electricity  occurs  via  absorption  of  light  by  the 
chromophore,  followed  by  the  interfacial  electron  transfer  (IET)  between  the  chromophore  and 
semiconductor.  The  most  successful  class  of  chromophores  used  in  the  molecular  assemblies  for  solar 
energy  conversion  is  based  on  Ru(II)-polypyridine  compounds.1  Unfortunately,  ruthenium  is  a  rare  and 
expensive  metal,  which  is  one  of  the  limiting  factors  to  its  wide  application  in  the  DSSCs  or  in 
photocatalytic  assemblies. 

Several  first  row  transition  metal  complexes  have  been  investigated  as  potential  photosensitizers, 
among  them  iron  and  copper  coordination  compounds.2  4  Our  studies  focus  on  Fe(II)  polypyridine 
complexes  due  to  their  close  relationship  with  the  very  successful  Ru(II)-polypyridines.  Additionally, 
Fe(II)  coordination  compounds  have  been  studied  in  their  role  as  photosensitizers  experimentally  by 
several  research  groups,3"8  which  provide  us  with  a  starting  point  and  a  set  of  benchmarks  for  our 
theoretical  calculations. 

In  this  report,  we  describe  the  results  of  our  efforts  to:  (1)  Evaluate  DFT  as  a  tool  for  calculations 
of  Fe(II)  complexes  and  show  that  it  is  an  adequate  tool  to  study  their  ground  and  excited  state  properties. 
(2)  Obtain  a  better  understanding  of  the  previous  experimental  findings  on  Fe(II)  complexes 
investigated  as  photosensitizers  and  published  in  the  literature.3"8  Our  results  provide  an  explanation  for  the 
band-selective  electron  injection  phenomena  in  Fe(II)-polypyridine-Ti02  nanoparticle  assemblies  and 
provide  a  pathway  toward  the  design  of  more  efficient  Fe(II)-based  chromophores. 
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3.  List  of  Appendices,  Illustrations,  and  Tables 


Table  1.  Experimental  and  calculated  spin  multiplicity  of  pseudo-octahedral  first  row  transition  metal 
complexes  investigated  in  this  work.  *  denotes  thermal  spin  crossover  complexes. 

Table  2.  Table  of  average  change  in  metal  -  ligand  bond  lengths  AR( metal- ligand)  in  Angstroms  for  low 
spin  to  high  spin  transition  versus  slope  of  the  scan  over  the  exact  exchange  for  various  pseudo-octahedral 
complexes.  *  [0s(bpy)3]2+  quintet  changes  electronic  states  over  the  exact  exchange  admixture 
investigated. 

Table  3.  Table  of  metal  ligand  bond  lengths  for  the  complexes  investigated,  coordination  site  labeling 
corresponds  with  Figure  6. 

Figure  1.  Fe(II)  complexes  investigated  in  this  work. 

Figure  2.  Heteropyridine  ligands  considered  in  this  study. 


Figure  3.  Dependence  of  the  energy  difference,  AEHS/LS  =  Ehj  h  in  -  Elow_spin  (kcal/mol),  on  the  exact 

exchange  admixture  (c,)  in  the  B3LYP  functional  for  iron  pseudo-octahedral  complexes  with  various 
ligands:  [Fe(CN)6]4-  (1),  [Fe(bpy)2(CN)2]°  (2),  [Fe(bpy)3]2+  (3),  [Fe(NH3)6]3+  (4),  [FeCl6f  (5).  A EHS/LS>  0 
corresponds  to  the  low-spin  ground  state  (singlet  or  doublet),  A EHS/ls<  0  corresponds  to  the  high-spin 
ground  state  (quintet  or  sextet). 


Figure  4.  Plot  of  average  change  in  metal  to  ligand  bond  lengths  in  transition  from  low-spin  to  high-spin 
states  versus  slope  of  the  scan  over  exact  exchange  for  data  shown  in  Table  2  along  with  a  plot  of  the 
linear  regression,  R2=0.93. 


Figure  5.  Dependence  of  the  energy  difference,  A EHSILS  =  Ehigh_spin-Elow_spinAEHS/LS  =  Ehigh_spin- 

Eiow-spin  (kcal/mol),  on  the  exact  exchange  admixture  (c/)  in  the  B3LYP  functional  for  several  pseudo- 
octahedral  iron(II)-polypyridine  complexes:  [Fe(bpy)2(CN)2]°  (1),  [Fe(bpy)3]2+  (2),  [Fe(tren(py)3)]2+  (3), 
[Fe(bpy)2(NCS)2]0  (4),  [Fe(tren(6-Me-py)3)]2+  (5),  and  [Fe(bpy)2Cl2]°  (6).  A Ehs/ls  >  0  A EHS/LS  >  0 
corresponds  to  the  singlet  ground  state,  A EHS/ls  <  0  AEHS/LS  <  0  corresponds  to  the  quintet  ground  state.  A 
linear  regression  is  plotted  for  each  complex  based  on  the  constant  slope  formula  A EHS/LS  =  I  -  140.2c, . 

Figure  6.  Optimized  structures  of  dyes  1,2,  and  3  with  each  coordination  site  labeled  (L1-L6). 

Figure  7.  Simulated  absorption  spectra  for  (from  top  to  bottom)  dyes  1,  2,  and  3  with  Lorentzian 
broadening  of  HWMH  =  0.12  eV  with  important  excitations  (fosc  >  0.01)  in  the  visible  region  labeled, 
calculated  with  B3LYP  TD-DFT  in  PCM  (acetonitrile). 

Figure  8.  Relevant  particle  states  for  the  major  excitations  (fosc  >  0.01)  from  the  calculated  spectra  shown 
in  Figure  7.  Kohn-Sham  orbitals  are  classified  by  the  absorption  band  (columns)  and  by  similarity  in  nodal 
structure  (rows). 

Figure  9.  Optimized  structures  of  anatase  (101)  with  pyridine-4-carboxylic  acid  binding  via  a 
monodentate  (a)  and  bidentate  (b)  binding  mode,  and  hydrogen  cyanide  (c). 
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Figure  10.  Final  slab  models  of  dye-nanoparticle  assemblies  consisting  of  [Fe(bpy-dca)2(CN)2]  attached  to 
anatase  (101)  via  a  monodentate  carboxylic  acid  binding  mode  in  the  following  nonequivalent 
orientations:  bpy  parallel  (left)  and  bpy  perpendicular  (right). 

Figure  11.  Structure  of  1-Ti02  assembly  employing  CN'  as  the  anchoring  group. 

Figure  12.  Density  of  states  for  1-Ti02  anatase  (101)  slab  model  (left)  and  enlarged  conduction  band 
(right)  show  the  following:  total  density  of  states  (blue  line),  projected  density  of  states  on  the  dye  (black 
line),  and  the  energy  levels  of  the  dye  in  vacuum  using  extended  Hiickel  (black  level  set  lines).  Gaussian 
line-shape  (HWHM  =  0.05  eV)  used  for  convolution. 

Figure  13.  Discrete  energy  levels  of  the  dyes  (left)  in  vacuum,  using  extended  Hiickel  level  of  theory, 
showing  the  evolution  of  the  similar  molecular  orbitals,  going  from  the  smallest  to  the  largest  dye.  The 
DOS/pDOS  for  the  [Fe(bpy-dca)2(CN)2]/Ti02  nanoparticle  system  is  shown  for  reference  (right). 

Figure  14.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of  [Fe(bpy)(CN)4]2' 
attached  to  anatase  (101)  slab  in  both  the  “bpy  parallel”  and  “bpy  perpendicular”  orientations  via 
monodentate  cyanide  anchoring  groups  with  varying  nitrogen  titanium  distances  (2.33,  2.22,  and  2.12  A), 
determined  by  either  exponential  or  biexponential  fitting. 

Figure  15.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of  [Fe(bpy- 
dca)(CN)4]2"  attached  to  anatase  (101)  slab  by  monodentate  carboxylic  acid,  bidentate  carboxylic  acid,  and 
cyanide  anchoring  groups  (in  both  the  “bpy  parallel”  and  “bpy  perpendicular”  orientations),  determined  by 
either  exponential  or  biexponential  fitting. 

Figure  16.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of 
[Fe(bpy-dca)2(CN)2]  attached  to  anatase  (101)  slab  by  monodentate  carboxylic  acid  (in  both  the  “bpy 
parallel”  and  “bpy  perpendicular”  orientations),  bidentate  carboxylic  acid  (in  both  the  “bpy  parallel”  and 
“bpy  perpendicular”  orientations),  and  cyanide  anchoring  groups,  determined  by  either  exponential  or 
biexponential  fitting. 
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4.  Statement  of  the  Problem  Studied 

Iron  pseudo-octahedral  complexes  have  been  long  considered  ideal  building  blocks  for  molecular 
electronic  switches,  data  storage  materials,  or  display  devices.9  Iron,  as  a  first  row  transition  metal, 
provides  a  weak  ligand  field  in  a  pseudo-octahedral  environment,  which  leads  to  the  presence  of  two 
energetically  close  electronic  states,  *A  and  5T.  Depending  on  the  character  of  ligands  coordinated  to  the 
central  iron  metal,  the  complexes  can  display  either  a  low-spin  ( 1  A),  or  a  high-spin  (5T)  ground  state.  The 
presence  of  a  low-lying  excited  state  of  a  different  spin  than  the  ground  state  is  also  responsible  for  the 
spin  crossover  phenomenon  in  these  compounds,  in  which  the  complex  changes  its  spin  state  under  the 
application  of  an  external  perturbation  such  as  a  change  in  temperature,  pressure,  or  exposure  to 
electromagnetic  radiation.10 

Fe(II)-polypyridine  complexes  have  also  been  investigated  as  potential  photosensitizers  in  dye- 
sensitized  solar  cells  due  to  their  structural  resemblance  to  the  Ru(II)-polypyridine  dyes  and  the  low  cost 
and  low  toxicity  of  iron.11-16  The  ability  of  Fe(II)  polypyridines  to  serve  as  photosensitizers  in  DSSCs  was 
first  demonstrated  by  Ferrere  and  Gregg,  although  their  efficiency  is  much  lower  than  that  of  their  Ru 
analogs.13  While  they  initially  absorb  visible  light  into  similar  metal-to-ligand  charge-transfer  (MLCT) 
states  as  Ru(II)-polypyridines,  the  photoactive  MLCT  states  undergo  intersystem  crossing  (ISC)  into  the 
low-lying  high-spin  metal-centered  states  on  sub-picosecond  time  scale.16  Therefore,  the  interfacial 
electron  transfer  (IET)  between  the  excited  dye  and  semiconductor  occurs  only  from  the  initially  populated 
“hot”  MLCT  states.  Since  the  lifetime  of  the  MLCT  manifold  of  these  systems  is  approximately  100  fs, 16-18 
the  IET  must  occur  with  a  characteristic  time  of  approximately  100  fs  or  less  in  order  to  be  competitive 
with  the  ultrafast  ISC  events. 

The  work  of  Ferrere  and  Gregg  also  established  the  band-selective  behavior  of  the  IET  between 
the  short-lived  excited  MLCT  states  of  Fe(bpy-dca)2(CN)2  (bpy-dca  =  2,2’-bipyridine-4,4,-dicarboxylic 
acid)  and  TiCL  semiconductor.13  Ferrere  et  al.  observed  that  upon  the  initial  excitation  with  visible  light, 
the  IET  is  more  efficient  (10-1 1%)  from  the  higher  energy  MLCT  band  and  much  less  efficient  (2%),  from 
the  lower  energy  MLCT  transitions. 

Computational  studies  of  Fe(II)-complexes  aiming  to  either  obtain  a  deeper  understanding  of 
their  electronic  structure  or  to  suggest  new  compounds  with  desirable  properties  face  a  number  of 
challenges.  One  of  these  challenges  is  the  correct  determination  of  the  ground  state,  which  is  non  trivial.19" 
20  There  have  been  a  number  of  computational  studies  dedicated  to  spin  crossover  compounds  using  the 
density  functional  theory  (DFT)20-28  as  well  as  higher  levels  of  theory,  such  as  CASPT2  (complete  active 
space  with  second-order  perturbation  theory).29  30  Due  to  the  size  of  these  systems  (50  or  more  atoms), 
they  are  most  amenable  to  calculations  performed  by  the  DFT  methods.  While  DFT  is  very  successful  at 
predicting  geometries  of  both  high-spin  and  low-spin  complexes,19  obtaining  the  correct  ground  state 
represents  a  major  challenge  since  GGA  (generalized  gradient  approximation)  functionals  (e.g.,  PBE)  tend 
to  favor  the  low  spin  states,  while  the  hybrid  functionals  (e.g.,  B3LYP)  artificially  favor  the  high  spin 
states.  Reiher  and  coworkers  related  this  behavior  to  the  amount  of  exact  Hartree-Fock  exchange  used  in 
the  hybrid  functionals.25  While  it  is  difficult  to  accurately  determine  the  ground  state  of  spin  crossover 
complexes  or  the  exact  energy  differences  between  their  high-spin  and  low-spin  states  with  the  DFT 
methods,  it  is  at  least  possible  to  qualitatively  predict  the  effect  of  ligand  substitution  on  the  spin  transition 
behavior.20’ 23-24 

The  objectives  of  this  work  are  to:  (1)  evaluate  DFT  as  a  tool  for  calculations  of  Fe(II)  complexes 
in  their  function  of  photosensitizers,  (2)  obtain  a  better  understanding  of  the  experimental  findings  on 
Fe(II)  complexes  investigated  as  photosensitizers  and  published  in  the  literature.3-8 

The  description  of  our  results  is  organized  into  two  parts  as  follows:  First,  we  summarize  our 
results  on  the  applicability  of  the  B3LYP  functional  to  the  ground  state  determination  of  pseudo- 
octahedral  iron  complexes,  with  special  focus  on  Fe(II)-polypyridines.  We  explore  systematic  trends  in  the 
calculated  energy  differences  between  the  high-spin  and  low-spin  electronic  states  for  a  series  of 
experimentally  known  transition  metal  complexes.  Finally,  we  suggest  that  the  ground  state  of  an  arbitrary 
Fe(II)  complexes  can  be  determined  by  comparing  the  calculated  energy  difference  between  the  low-spin 
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and  high-spin  electronic  states  to  the  energy  differences  of  structurally  related  complexes  with  a  known, 
experimentally-determined  ground  state,  without  the  need  to  reparametrize  the  functional. 

Second,  we  describe  our  studies  on  the  light  absorption  and  IET  processes  for  three  different  dye- 
TiC>2  nanoparticle  systems,  with  a  special  focus  on  elucidating  the  origin  of  the  band-selective  IET  in 
Fe(II)  polypyridine-TiCE  assemblies.  The  dyes  investigated  include  [Fe(bpy)(CN)4]2'  (1),  [Fe(bpy- 
dca)(CN)4]2'  (2),  and  [Fe(bpy-dca)2(CN)2]  (3),  shown  in  Figure  1.  Two  of  these  dyes  were  studied 
previously,  1  by  Meyer14'15  and  3  Ferrere.11' 13  While  3  attaches  to  the  Ti02  surface  via  the  carboxylic  acid 
anchoring  group,  1  is  attached  via  the  CN'-Ti  bridge.  Although  dye  2  was  not  studied  experimentally,  it 
represents  a  natural  bridge  between  the  complexes  1  and  3  and  was  therefore  included  in  our 
computational  study. 

5.  Summary  of  the  most  important  results 


5.1.  Benchmarking  the  DFT  methodology 
5.1.1  Computational  Approach 

Twenty  seven  pseudo-octahedral  complexes  of  Fe,  Co,  Mn,  Ru  and  Os  with  a  mixture  of 
heteropyridine  and  Cl',  F',  CN‘,  NCS",  and  H20  ligands  were  investigated.  The  majority  of  the  complexes 
(19)  are  complexes  of  iron.  Heteropyridine  ligands  considered  are  shown  in  Figure  2.  They  include  bpy  = 
2,2'-bipyridine,  phen  =  1 , 1 0-phenanthroline,  terpy  =  2,2';6',2"-terpyridine,  tren(py)3  =  tris((N-(2- 
pyridylmethyl)-2-iminoethyl)amine),  tren(6-Me-py)3  =  tris((N-(2-(6-Me-pyridyl)methyl)-2- 
iminoethyl)amine),  Lt  =  4-(4-bromophenyl)-2,6-bis(2-pyridine)-l,3,5-triazine,  Me-Lt  =  4-(4- 

bromophenyl)-2,6-bis(6-picolinyl)-l,3,5-triazine,  tpen  =  tetrakis(2-pyridylmethyl)ethylenediamine,  and 
(pyrol)3tren  =  tris(l-(2-azoyl)-2-azabuten-4-yl)amine  ligands. 

Geometries  of  all  compounds  were  initially  optimized  at  the  B3LYP  level  of  theory.31  SDD 
relativistic  effective  core  potential  and  associated  triple-^  basis  set32  was  used  to  describe  the  metal  center 
and  6-31 1+G*  basis  set33'34  was  used  for  all  other  atoms.  Singlet  and  quintet  states  were  optimized  for 
complexes  of  Fe(II),  Co(III),  Ru(II),  and  Os(II),  doublet  and  sextet  states  were  optimized  for  complexes  of 
Fe(III)  and  Co(II),  triplet  and  quintet  states  were  optimized  for  the  Mn(III)  complex. 

To  investigate  the  dependence  of  high-spin  vs.  low-spin  energy  differences 

(  HS/LS  =  E high-spin  ~  ^ low-spin  >  HS/LS  =  H high-spin  ~  ^ low-spin  ’  HS/LS  =  ^ high-spin  ~  ^ low-spin  )  011  l^*e  ^mOUnt  of 


exact  exchange  in  the  B3LYP  functional,  we  have  systematically  varied  c/  parameter  from  0.0  to  0.25, 
corresponding  to  0-25%  Hartree-Fock  exchange  in  the  exchange-correlation  functional  of  B3LYP: 


c3AE£W91,  (1) 

in  which  c/,  c2,  and  c.;  are  coefficients  fit  from  experimental  data  (determined  to  be  c/=  0.20,  c2  =  0.72,  c3 
=  0.81),31  Exlcsda  is  the  exchange-correlation  energy  from  the  local  spin  density  approximation  (LSD A), 35 
Eexact  is  the  exact  exchange  energy,  2?f88  is  Becke’s  1988  gradient  correction  of  exchange,36  and  A E™91 

is  Perdew  and  Wang’s  1991  gradient  correction  to  correlation.37 

The  low  spin  and  high  spin  geometries  of  all  compounds  were  fully  re-optimized  for  each  different 
value  of  the  exact  exchange  in  the  B3LYP  functional.  All  optimized  structures  were  verified  using 
vibrational  frequency  analysis.  The  results  of  vibrational  analysis  were  used  to  obtain  enthalpies  ( H )  and 
free  energies  ( G )  for  all  compounds  at  298.15  K.  An  ultra- fine  integration  grid  was  used  for  all 
calculations. 

Calculated  expectation  values  of  the  S2  operator  are  within  10%  of  the  value  expected  for  each 
complex,  except  for  the  doublet  spin-state  of  [FeF6]3'  (<  S2>  =  1.79  in  vacuum  and  <  S2>  =  1.14  in  water) 
and  [FeCl(,]3'  (<  S2>  =  1.67  in  vacuum  and  <  S2>  =  1.30  in  water).  The  resulting  wavefunctions  were  also 
tested  for  stability.38'39  All  species  were  optimized  in  vacuum  and  solvent  (water)  using  Polarizable 
Continuum  Model  (PCM),  with  the  exception  of  [Fe(Me-Lt)2]2+  which  was  relaxed  in  vacuum  only.40  The 
use  of  solvent  improved  the  structures  and  spin  contamination  of  anionic  compounds,  while  it  did  not  have 
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any  significant  influence  on  optimized  geometries  of  cations  and  neutral  species.  All  calculations  were 
performed  using  Gaussian  09  software  package.41 

5.1.2.  Results  and  Discussion 

Ground  state  determination  by  B3LYP,  energies  of  high-spin  vs.  low-spin  states.  While 
B3LYP  can  be  used  as  a  reliable  tool  in  obtaining  high-spin  and  low-spin  geometries  of  first  row  transition 
metal  complexes,  it  is  not  as  reliable  in  predicting  their  ground  state  spin  multiplicity.  Out  of  the  25  first 
row  transition  metal  complexes  investigated,  B3LYP  was  incorrect  in  predicting  the  ground  state  of  eight 
of  these,  artificially  favoring  the  high-spin  state  over  the  low-spin  state  in  seven  of  the  cases  (see  Table  1). 
Inclusion  of  thermal  and  entropic  factors  further  stabilized  high-spin  state  in  four  more  complexes  and 
incorrectly  predicted  the  high-spin  ground  state  for  10  of  the  25  complexes.  Such  bias  toward  the  high-spin 
states  is  not  unusual  for  the  B3LYP  and  has  been  observed  before.26 

Following  the  example  of  Reiher  et  al., 26  we  have  calculated  the  energy  (A EHS/ls),  as  well  as 
enthalpy  (A HHS/ls)  and  free  energy  differences  (A GHs/ls)  between  the  high-spin  and  low-spin  states  for  a 
group  of  transition  metal  compounds  with  respect  to  the  exact  exchange  admixture  in  the  B3LYP 
functional  form  (c/  =  0.00  -  0.25,  see  Equation  1).  The  energy  difference  was  defined  as 

AE HS/LS  =  E high-spin  ~  E low-spin  (  AH HS/LS  ~  H high-spin  ~  H low-spin  ’  AG HS/LS  ~  Er high-spin  ~  ^ low-spin  )  &nd  POSitive  if 

the  low-spin  state  is  calculated  to  be  more  stable  than  the  high-spin  state. 

Figure  3  shows  the  AEhs/ls  dependence  on  a  for  a  number  of  different  Fe(II)  pseudo-octahedral 
complexes  with  polypyridine  as  well  as  non-polypyridine  ligands.  While  AEHS/ls  depends  linearly  on  c/, 
the  slope  of  this  dependence  varies  widely  among  the  complexes  and  it  is  not  possible  to  fit  this 
dependence  with  a  single  linear  equation  for  all  complexes  investigated.  Similar  behavior  has  been 
reported  previously  by  Salomon  et  al.27  for  a  different  set  of  first  row  transition  metal  complexes. 
Interestingly,  [FeF6]3'  displays  nonlinear  behavior  in  its  dependence  of  A EHS/ls  on  c/  when  optimizations 
are  performed  in  vacuum.  Further  investigation  of  this  phenomenon  revealed  the  presence  of  two  different 
high-spin  states  for  this  complex,  whose  relative  stabilities  change  depending  on  the  amount  of  exact 
exchange  in  the  B3LYP  functional.  Therefore,  calculations  performed  with  lower  admixture  of  exact 
exchange  converge  to  a  different  electronic  state  than  those  with  a  higher  admixture,  resulting  in  the 
nonlinear  behavior.  This  problem  of  state  switching  for  [Fe(F)6]3'  is  not  present  when  optimizations  are 
performed  in  solvent  (FEO)  using  PCM. 

Table  2  summarizes  the  behavior  of  all  complexes  investigated,  giving  the  slope  of  the  linear  fit 
for  the  dependence  of  AEHS/ls,AHHS/LS,  and  AGHs/ls  on  the  amount  of  exact  exchange  {ci).  Interestingly,  all 
Fe(II)-polypyridine  complexes  have  similar  slope  of  A E  vs.  c/  across  the  0.0  -  0.25  range  of  the  exact 
exchange  admixture.  On  the  other  hand,  while  the  dependence  of  A E  vs.  c/  for  more  ionic  complexes  (e.g., 
[Fe(H20)]2+,  [FeCh,]4")  is  still  linear,  the  slope  of  this  dependence  is  completely  different  from  that  of  the 
Fe(II)-polypyridine  compounds.  This  makes  it  challenging  to  come  up  with  a  single  correction  scheme  for 
a  wide  variety  of  transition  metal  compounds.  We  suggest  that  there  is  a  connection  between  the  slope  of 
the  A E  vs.  ci  dependence  and  structural  features  of  the  complexes,  as  well  as  the  character  of  metal-ligand 
bonds  in  the  complex  (ionic  vs.  covalent). 

Structure  of  Fe(II)  complexes  and  error  in  high-spin/low-spin  energies  from  B3LYP. 

According  to  the  ligand  field  theory,  the  ligand  field  splitting  between  a  set  of  t2g  and  eg  metal  orbitals  is 
determined  by  the  identity  of  a  central  atom  and  its  ligands.  The  size  of  the  ligand  field  splitting 
determines  the  ground  state  electronic  configuration  (high-spin  vs.  low-spin)  as  well  as  the  energy 
difference  between  the  high-spin/low-spin  states  of  a  particular  octahedral  compound.  Therefore,  it  would 
not  be  surprising  if  there  were  some  link  between  the  structure  of  the  complexes  investigated  and  the  slope 
of  AEhs/ls  vs.  ci. 

The  change  in  the  electronic  configuration  from  high-spin  to  low-spin  and  vice  versa  is  usually 
accompanied  by  a  marked  change  in  the  metal-ligand  bond  lengths  of  the  covalent  complexes,  as  placing 
the  electrons  into  the  antibonding  eg  orbitals  increases  the  metal-ligand  distance  in  the  high-spin  state. 
Therefore,  we  have  explored  the  relationship  between  the  average  change  in  the  metal-ligand  distance 
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(AR)  and  the  slope  of  the  A EHs/ls  dependence  on  amount  of  exact  exchange  in  the  B3LYP  functional  (<?/). 
The  reported  AR  is  an  average  over  the  ci  range  investigated,  but  the  results  are  virtually  identical  when 
only  considering  AR  calculated  at  cj  =  0.2  (unmodified  B3LYP  functional). 

Table  2  summarizes  the  average  change  in  the  metal-ligand  distance  between  the  high-spin  and 
low-spin  state  determined  at  the  B3LYP  level  of  theory  along  with  the  slope  of  the  A Ehs/ls,  AHHS/ls,  and 
A Ghs/ls  dependence  on  cj.  The  change  in  the  metal-ligand  distance  was  determined  as  an  average  of  the 
changes  calculated  for  all  values  of  c \  investigated  between  0.0  and  0.25.  Figure  4  shows  the  plot  of  the 
average  change  in  the  metal-ligand  bond  lengths  between  the  high-spin  and  low-spin  states,  with  respect  to 
the  slope  of  the  A Ehs/ls  vs.  cj  dependence.  Interestingly,  there  is  a  linear  relationship  between  the  two, 
showing  that  a  smaller  slope  corresponds  to  a  smaller  change  in  the  bond  lengths.  The  R2  coefficient  for 
this  dependence  is  0.93,  suggesting  a  strong  correlation.  Second  and  third-row  transition  metal  compounds 
([Ru(bpy)3]2+,  [0s(bpy)3]2+,  shown  in  blue)  were  excluded  from  the  fit,  but  are  still  shown  on  the  plot.  A 
weaker  correlation  is  observed,  with  an  R2  value  of  0.78,  when  all  species  are  optimized  in  vacuum.  All 
species  reported  in  Table  2  were  optimized  using  PCM  solvent  model  for  H20  excluding  [Fe(Me-Lt)2]2+, 
which  was  optimized  in  vacuum. 

It  is  apparent  that  the  complexes  with  the  weakest  dependence  of  A Ehs/ls  on  c/  are  ionic 
compounds  of  first  row  transition  metals  (e.g.,  [Fe(H20)6]2+,  [FeCf,]  ’")  that  also  undergo  smallest  change 
in  the  metal-ligand  bond  lengths  between  the  low  spin  and  high  spin  states.  Complexes  with  more  covalent 
character  of  metal-ligand  bonds  undergo  larger  structural  changes  between  their  low-spin  and  high-spin 
electronic  states  and  are  more  strongly  influences  by  the  amount  of  exact  exchange  admixture  in  the  DFT 
functional. 

The  outliers  to  this  trend  are  shown  in  blue  and  belong  to  octahedral  compounds  with  a  central 
atom  from  second  and  third  row  of  transition  metals  (Ru  and  Os),  suggesting  that  the  relationship  between 
the  average  change  in  metal-ligand  bond  lengths  and  the  slope  of  the  A Ehs/Ls  dependence  on  c2  will  have  a 
different  character  than  displayed  by  the  complexes  of  first  row  transition  metals. 

Note  that  a  similar  relationship  as  observed  between  the  dependence  of  A EHS/LS  on  c/  vs.  AR  holds 
true  for  AHhs/Ls  and  AGhs/ls ■  However,  the  correlation  for  this  dependence  becomes  slightly  worse  with  the 
thermal  and  entropic  corrections  included  (R2=  0.93  and  0.91,  respectively).  This  is  not  very  surprising,  as 
several  assumptions  are  made  in  obtaining  these  corrections  (i.e.,  assumption  of  harmonic  potential, 
neglecting  the  role  of  the  solvent)  that  are  not  equally  good  for  all  compounds  investigated  in  this  study. 
The  problem  is  exacerbated  by  the  fact  that  the  low  frequency  vibrations,  which  provide  the  largest 
contribution  to  the  entropy,  suffer  from  the  highest  error  in  the  harmonic  approximation.25 

It  is  also  worth  mentioning  that  identical  results  are  obtained  when  the  analysis  is  performed  with 
a  smaller  basis  set  (6-31G*  instead  of  6-31 1+G*)  and  in  vacuum.  The  larger  basis  set  and  solvent  model 
are,  however,  necessary  in  order  to  obtain  proper  description  of  the  structures  and  electronic  states  of  some 
anionic  compounds  (e.g.  [Fe(CN)6]4',  [FeF6]3')- 

Ground  state  determination  for  a  family  of  structurally  related  complexes.  The  linear 
correlation  between  the  average  AR  and  slope  of  the  A EHS/ls  vs.  c/,  suggests  that  the  octahedral  complexes 
that  undergo  similar  distortion  in  their  octahedral  environment  going  from  a  low-spin  to  high-spin  state 
(described  here  by  the  average  change  in  metal-ligand  bond  lengths),  suffer  from  the  similar  error  in  the 
B3LYP  functional.  In  case  of  the  covalent  complexes  of  iron,  this  systematic  error  results  in  the  artificial 
stabilization  of  their  high-spin  state  energies  with  respect  to  the  low-spin  state  energies.  The  question  still 
remains  if  this  error  is  systematic  for  a  group  of  similar  complexes,  i.e.  if  it  results  in  the  artificial 
stabilization  of  the  high-spin  state  energies  with  respect  to  the  low-spin  state  energies  by  about  the  same 
amount. 

To  find  if  the  error  is  indeed  systematic,  we  have  looked  at  the  A EHS/LS  vs.  c2  dependence  in  a  set  of 
ten  structurally-related  Fe(II)-polypyridine  complexes.  They  were  chosen  as  “structurally-related”  because 
they  display  similar  geometry  change  in  the  metal-ligand  bond  lengths  going  from  low  spin  to  high  spin 
geometry  (see  Table  2).  The  selected  complexes  contain  polypyridine  ligands  combined  with  other  ligands 
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of  varying  ligand-field  strength,  such  as  Cl",  NCS",  and  CN".  The  plot  of  the  AEks/ls  as  a  function  of  ci  for 
six  members  from  this  family  of  related  Fe(II)  compounds  is  shown  in  Figure  5. 

As  can  be  seen  in  Figure  5,  the  dependence  of  A EHS/LS  on  a  is  linear,  with  approximately  the  same 
slope  for  the  compounds  considered.  Note  that  very  similar  linear  dependence  of  A EHs/ls  on  a  as  well  as 
systematic  errors  for  a  group  of  related  Fe(II)-S  complexes  have  been  observed  previously  by  Reiher  and 
coworkers.26  In  case  of  our  Fe(II)  test  complexes,  this  linear  relationship  between  A E  and  c/  can  be 
described  by  the  following  equation: 

Mhsils  =  /-  142.5c,  AEhs/ls  =  I-140.2cls  (2) 

where  /  corresponds  to  the  intercept  (the  value  of  A Ehs/ls  at  cj=  0  for  each  complex)  which  describes  the 
vertical  shift  between  the  plots  of  A EHS/ls  vs.  cj  for  different  complexes.  The  four  complexes  included  in 
the  fit  but  not  shown  in  Figure  5  are:  [Fe(Lt)2]2+,  [Fe(Me-Lt)2]2+,  [Fe(phen)3]2+,  [Fe(terpy)2]2+,  and 
[Fe(tpen)]2+.  The  above  equation  fits  the  AEhs/ls  dependence  on  a  for  each  of  the  eleven  complexes  very 
accurately,  with  the  correlation  coefficient  R1  =  0.97-0.99.  More  importantly,  the  vertical  shift  between 
plots  of  AEHs/Ls  for  different  Fe(II)  compounds  correctly  reflects  the  change  in  the  ligand  field  strength  of 
their  ligands,  meaning  that  the  low-spin  states  of  Fe(II)  compounds  are  stabilized  over  the  high-spin  states 
in  the  same  order  as  the  Fe(II)  ligands  appear  in  the  spectrochemical  series:  Cl"  <  NCS"<  pyridine  < 
bipyridine  <  CN".  This  means  that  although  B3LYP  is  not  capable  to  predict  the  correct  ground  state  spin 
multiplicities,  the  overall  order  of  the  high-spin/low-spin  energy  differences  for  a  group  of  structurally 
related  compounds  reflects  their  high-spin,  low-spin,  or  thermal  spin  crossover  character  in  the  ground 
state. 

It  should  be  noted  that  when  comparing  [Fe(phen)3]2+  and  [Feibpyf]2  we  found  that  the  trend  in 
AEhs/ls  vs.  c2  incorrectly  reflected  the  trend  in  ligand  field  strength,  however  the  distance  between  the  two 
plots  of  A Ehs/ls  over  c/  =  0-0.25  were  at  most  2.0  kcal/mol.  [Fe(phen)3]2+ is  not  shown  in  Figure  5,  but 
belongs  to  the  same  set  of  structurally-related  complexes.  Phen  and  bpy  ligands  are  considered  to  have 
very  similar  ligand  field  strength,  so  this  discrepancy  is  well  within  the  accepted  error  range  for  the 
B3LYP  functional. 

The  linear  behavior  of  A EHS/ls  with  respect  to  c/,  along  with  the  systematic  trend  in  the  intercept  / 
for  a  series  of  related  compounds,  can  therefore  be  used  to  construct  a  benchmark  for  determination  of  the 
ground  state  multiplicity  for  an  arbitrary  Fe(II)-polypyridine  complex.  A  procedure  to  determine  the 
ground  state  multiplicity  of  an  arbitrary  Fe(II)-polypyridine  complex  could  be  as  follows:  (1)  Obtain 
A EHs/ls  vs.  exact  exchange  dependence  plots  for  a  series  of  related  compounds  with  available  experimental 
reference  data.  (2)  Based  on  this  benchmark,  decide  which  values  of  the  intercept  /  correspond  to  the  high- 
spin  and  low-spin  ground  states.  In  the  specific  example  of  Fe(II)-polypiridines  shown  in  Figure  5,  />  24 
kcal/mol  indicates  a  low-spin  complex,  /  <  18  kcal/mol  high-spin  complex,  with  /  =  18-24  kcal/mol 
indicating  a  possible  thermal  spin  crossover  complex.  (3)  Obtain  AEhs/ls  vs.  exact  exchange  dependence 
plot  for  the  compound  with  unknown  ground  state.  (4)  Based  on  the  intercept  I  obtained  for  the  unknown 
compound,  determine  its  spin  state. 

For  a  series  of  related  complexes  this  procedure  will  produce  a  unique  range  of  /  values  for  low- 
spin  and  high-spin  stability.  A  questionable  spin  crossover  region  will  also  always  exist  between  the  two 
stability  regions  due  to  inherent  error  in  the  method.  The  value  of  /  for  a  thermal  spin  crossover  complex 
determines  where  the  low-spin  and  high-spin  states  change  relative  stability.  The  average  of  the  intercepts 
for  the  two  thermal  spin  crossover  complexes  ([Fe(tpen)]2+,  [Fe(bpy)2(NCS)2]°)  is  21  kcal/mol  which  is, 
interestingly,  close  to  the  determinded  exchange  interaction  between  electrons  in  the  d  orbitals  of  iron.42"44 
We  believe  this  is  a  coincidence  since  the  intercept  for  a  spin  crossover  complex  represents  the  error  in 
AEHs/ls  when  the  B3LYP  functional  is  reparameterized  to  be  a  purely  GGA  functional  (ci=0),  which  still 
contains  exchange  interaction. 

The  strong  linear  correlation  between  the  average  A R  and  slope  of  the  AEhs/ls  vs.  c/  scans  could 
also  allow  us  to  forgo  the  calculation  of  A EHs/ls  vs.  c2  dependence  for  several  values  of  c2  for  both 
benchmark  and  unknown  complexes.  It  makes  it  possible  for  us  to  rely  on  the  average  change  in  the 
octahedral  coordination  environment  to  determine  what  group  of  complexes  is  “similar”  for  the  purposes 
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of  the  ground  state  determination  and  reduces  the  number  of  calculations  one  needs  to  perform  to  just 
B3LYP  with  a  single  value  of  c;  =  0.2.  Analogous  analysis  with  virtually  identical  results  toward  lower 
energies  can  also  be  performed  using  AHhs/Ls  or  A GHs/ls  instead  of  A EHS/ls,  noting  that  the  calculated 
A Hhs/ls  and  AG/m,s  will  be  systematically  shifted  toward  lower  energies. 

5.2.  Elucidating  band-selective  injection 
5.2.1  Computational  Approach 

Molecular  Structure  and  Absorption  Spectra.  Dye  molecules  were  optimized  at  the  B3LYP31'45 
level  of  theory.  The  SDD  relativistic  effective  core  potential  (ECP)  and  associated  basis  set32  were  used  to 
describe  the  central  iron  atom,  and  6-3 1G*  basis  sets33’34  were  used  for  all  other  atoms.  All  complexes 
were  optimized  in  the  singlet  ground  state  in  vacuum.  The  calculation  of  UV-vis  absorption  spectra  for 
dyes  1-3  was  performed  employing  TD-DFT46  48  methodology  with  the  same  basis  set  and  functional  as 
used  for  the  geometry  optimizations.  Only  spin  allowed  singlet  vertical  excitations  were  determined  and 
TD-DFT  was  performed  on  the  entirely  relaxed  dye  molecules  1-3.  Polarizable  continuum  model  (PCM)49, 
using  acetonitrile  as  a  solvent,  was  employed  in  the  TD-DFT  calculations.  Absorption  spectra  were 
simulated  by  convoluting  the  stick  spectrum  composed  of  the  8-functions  associated  with  each  excitation 
energy  times  the  oscillator  strength  with  a  Lorentzian  line  shape  with  half-width-at-half-maximum 
(HWHM)  of  0.12  eV.  Canonical  Kohn-Sham  orbitals  were  used  to  characterize  absorption  peaks  with  fosc 
>  0.01,  within  the  visible  portion  of  the  spectra  (A,  >  350  nm).  The  Gaussian  09  software  package41  was 
used  for  all  the  DFT  and  TD-DFT  calculations  on  dye  molecules. 

Slab  Model  Optimization.  The  Vienna  Ab-initio  Simulation  Package  (VASP)50  5j  was  used  to 
optimize  periodic  systems  (bulk  TiCF  and  nanoparticles)  at  DFT  level  of  theory.  The  Perdew-Burke- 
Ernzerhof  (PBE)54'55  exchange-correlation  functional  with  Projector  Augmented- Wave  method56  57  was 
employed  in  all  periodic  boundary  calculations  (PBCs).  The  plane  wave  basis  set  expansion  was  cutoff  at 
500  eV  for  all  PBCs.  The  unit  cell  for  anatase  TiCF  was  obtained  with  a  (13x13x1 3)  A-point  sampling, 
resulting  in  a  tetragonal  lattice  with  lattice  vectors  a  =  b  =  3.81  A,  c  =  9.77  A).  The  optimized  geometry  of 
the  TiC>2  unit  cell  was  used  to  construct  a  slab  model  of  the  (101)  TiC>2  surface  functionalized  with  either 
hydrogen  cyanide  or  pyridine-4-carboxylic  acid.  The  linker  models  and  the  top  2  layers  of  Ti  and  top  4 
layers  of  O  were  relaxed  using  a  (5  x  3  x  1)  /c-point  sampling,  keeping  the  supercell  volume  fixed  with 
lattice  vectors  a  =  15.25  A,  b  =  10.49  A,  c  =  26.00  A.  The  bottom  layers  of  the  slab  were  fixed  at  bulk 
geometry. 

The  final  model  dye-nanoparticle  assemblies  were  constructed  by  performing  a  constrained 
optimization  of  the  dyes  1-3  with  the  linker  groups  (cyanide,  and  carboxylic  acid  in  both  monodentate  and 
bidentate  surface  binding  modes)  fixed  at  surface-optimized  geometries.  Coordinates  for  the  carboxylate 
functional  group  as  well  as  the  three  closest  carbons  in  the  pyridine  ring  along  with  the  two  adjacent 
hydrogens  were  held  frozen  in  the  constrained  optimization  of  dyes  2  and  3.  The  Gaussian  09  software 
package,41  employing  the  same  functional  and  basis  sets  as  above,  was  used  for  these  constrained 
optimizations  of  dye  with  frozen  linker  geometry.  Dyes  2  and  3  were  attached  to  the  slab  surface  in  a 
cyanide,  monodentate  carboxylic  acid,  and  bidentate  carboxylate  binding  modes  via  the  Kabsch 
alignment.58  Dye  1  was  attached  in  the  same  manner  with  a  cyanide  binding  mode.  All  the  nonequivalent 
linker  positions  on  the  dyes  were  modeled:  two  nonequivalent  carboxylic  acid  groups  on  1  and  2  and  two 
nonequivalent  cyanide  groups  on  3.  All  optimizations  were  performed  in  vacuum. 

Interfacial  Electron  Transfer  Simulations.  All  model  systems  were  composed  of  a  (101)  anatase 
nanoparticle  slab  with  dyes  (1-3)  attached  via  an  anchoring  group  (cyanide  or  carboxylic  acid)  in  vacuum 
and  their  construction  is  described  above.  Quantum  dynamics  simulations  to  model  the  interfacial  electron 
transfer  (IET)  were  performed  on  each  dye-semiconductor  model  using  a  method  developed  by  Rego  and 
Batista,59  employing  the  extended  Hiickel  (EH)  Hamiltonian.  Details  of  this  method  can  be  found  in 
references  59  and  60.  Survival  probability  was  obtained  at  each  simulation  step  by  projecting  of  the  time- 
evolved  wave  function  onto  the  atomic  orbitals  of  the  adsorbed  chromophore.  Characteristic  electron 
injection  time  was  then  obtained  from  an  exponential  fit  of  the  survival  probability. 
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For  each  system,  donor  states  in  the  IET  simulations  were  chosen  from  the  particle  states  that 
make  up  major  contributions  to  prominent  excitations  in  the  visible  region  of  the  absorption  spectrum 
(fosc  >  0.01,  A  >  350  nm)  as  calculated  by  TD-DFT  calculations.  The  canonical  Kohn-Sham  orbitals 
corresponding  to  the  relevant  particle  states  were  matched  to  the  EH  orbitals. 

Simulations  were  run  with  a  time  step  of  0.1  fs  up  to  3000  fs  on  super-cell  arrays  with  30.49  A  x 
31.46  A  x  26.00  A  dimensions.  One  L-point  was  used  for  sampling  and  the  calculation  was  run  using 
periodic  boundary  conditions.  All  simulations  were  done  at  the  frozen  geometry.  Absorbing  potentials 
(i.e.,  imaginary  terms  to  the  diagonal  elements  of  the  Hamiltonian)  were  placed  on  the  bottom  layer  of  Ti 
atoms  to  avoid  artificial  recurrences  in  the  electron-transient  populations. 

The  survival  probability  curves  were  fit  with  either  a  exponential  decay  function,  or  in  the  case 
that  R2<  0.95,  a  biexponential  decay  function.  The  survival  probability  was  constrained  to  unity  at  time  0 
for  all  fits.  The  characteristic  injection  time  was  determined  by  the  reciprocal  of  the  decay  constant  for 
single  exponential  functions,  or  the  reciprocal  of  the  decay  constant  responsible  for  the  majority  of  decay 
(coefficient  >0.50)  for  the  biexponential  functions. 

5.2.2.  Results  and  Discussion 

The  results  are  organized  as  follows:  First  the  ground  state  molecular  properties  and  absorption 
spectra  are  described  for  all  dyes,  followed  by  the  description  of  the  dye-nanoparticle  assemblies.  Finally, 
results  of  the  IET  simulation  are  summarized  for  each  dye  individually. 

Ground  State  Structures  and  Absorption  Spectra.  The  optimized  structures  of  dyes  1  -  3  are 
shown  in  Figure  6.  Metal-ligand  bond  lengths  for  all  three  complexes  are  reported  in  Table  3,  using  the 
atom-numbering  scheme  from  Figure  6.  For  dye  1,  we  find  a  2%  error  when  comparing  B3LYP  metal- 
ligand  bond  length  to  the  experimental  crystal  structure  data.61  For  dye  2,  crystal  structure  data  is 
unavailable  but  we  assume  errors  are  similar  to  the  related  dye  1.  For  dye  3,  crystal  structure  data  is  only 
available  for  the  related  species  [Fe(bpy)2(CN)2],  When  comparing  the  metal-ligand  bond  lengths  for  the 
singlet  state  of  dye  3  to  the  corresponding  bonds  from  the  crystal  structure  of  [Fe(bpy)2(CN)2],  we  find 
only  a  2%  error.62  Overall,  B3LYP  provides  accurate  geometries  for  all  the  dyes. 

The  simulated  absorption  spectra  for  dyes  1-3  are  shown  in  Figure  7.  All  three  complexes  display 
two  absorption  bands  in  the  visible  region.  Excitations  A-C  (complex  1)  and  A-F  (complex  2)  can  be  all 
described  as  metal-to-ligand  charge  transfer  (MLCT)  transitions.  While  majority  of  the  transitions  A-J 
(complex  3)  also  have  a  MLCT  character,  several  transitions  (D,  E)  have  significant  metal-centered  (MC) 
character. 

Relevant  particle  states  involved  in  the  visible  light  transitions  are  shown  in  Figure  8,  classified  by 
the  absorption  band.  Note  that  functionalization  of  the  bipyridine  ligand  of  dye  1  with  dicarboxylic  acid, 
resulting  in  dye  2,  changes  the  energetic  ordering  of  the  LUMO+1  and  LUMO+2  orbitals.  Complex  3 
possesses  three  degenerate  sets  of  virtual  orbitals  (LUMO,  LUMO+1  and  LUMO+2,  LUMO+3)  and  a 
contribution  from  metal-centered  ligand  field  state  (LUMO+12).  The  particle  states  shown  in  Figure  8 
serve  as  initial  (donor)  states  for  the  IET  simulations  reported  below. 

Dye-Nanoparticle  Assemblies.  The  linker  models  used  to  attach  the  dye  molecules  to  the  Ti02 
surface  are  shown  in  Figure  9.  We  have  explored  cyanide  attachment  mode,  and  monodentate  and 
bidentate  attachment  modes  for  carboxylic  acid.  Although  multiple  carboxylic  acid  binding  modes  exist, 
only  two  attachment  modes  are  explored  in  this  work,  based  on  previously  published  stability 
analysis.63 

Each  of  the  Fe(II)  polypyridine  compounds  investigated  can  attach  to  the  surface  in  multiple  ways. 
Complex  3  can  be  anchored  to  the  Ti02  surface  via  -COOH  group  attached  at  4  or  4’  site  of  the  bpy 
ligand.  Similarly,  two  non- equivalent  CN'  sites  for  attachment  are  present  in  1  and  2.  This  results  in  two 
possible  attachments  of  all  three  dyes  onto  the  surface,  in  which  the  bipyridine  group  not  bound  to  the 
surface  is  oriented  either  parallel  or  perpendicular  relative  to  the  (101)  plane  of  Ti02  (see  Figure  10). 

Dyes  anchored  to  the  surface  via  the  CN'  group  have  a  free  degree  of  rotation  about  the  axis 
defined  by  the  cyanide  bond  (see  Figure  11).  This  ambiguity  is  caused  by  the  fact  that  HCN-Ti02  surface 
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model  was  chosen  to  model  the  full  dye-TiC>2  assembly  due  to  the  computational  constrains.  A 
conformation  that  places  the  bipyridine  groups  as  far  from  the  surface  as  possible  was  chosen  for  the  IET 
simulations  reported  here.  This  allows  us  to  study  the  IET  via  the  surface-anchoring  group,  rather  than 
indirect  IET  between  the  bpy  ligand  group  and  TiC>2  surface  due  to  their  spatial  proximity. 

It  is  important  to  note  that  the  structure  optimization  of  the  HCN/TiC>2  system  (see  Figure  11) 
resulted  in  a  significantly  longer  N-Ti  bond  length  (2.33  A)  than  that  reported  by  Selloni64  (2.12  A). 
Selloni  and  coworkers  have  optimized  [Fe(CN)6]4  on  a  Ti3g076  cluster,  which  likely  provides  a  better 
model  for  this  attachment.  The  influence  of  the  N-Ti  distance  on  the  results  on  the  IET  rates  was  therefore 
investigated  by  setting  the  N-Ti  bond  length  to  2.12,  2.22,  and  2.33  A  in  the  IET  simulations  as  shown  in 
Figure  11. 

Interfacial  Electron  Transfer  Simulations.  Density  of  states  plots  for  chromophore-nanoparticle 
assemblies  calculated  with  the  extended  Hiickel  method  are  shown  in  Figures  12  and  13.  For  all  three 
complexes  investigated,  dyes  introduce  occupied  energy  levels  into  the  semiconductor  band  gap  and  a 
number  of  virtual  energy  levels  (LUMO  -  LUMO  +  12  for  dye  1,  LUMO  -  LUMO  +  13  for  dye  2,  and 
LUMO  -  LUMO  +  15  for  dye  3)  that  couple  with  the  conduction  band  of  Ti02. 

Differences  in  the  orbital  energy  levels  among  the  three  dyes  are  illustrated  in  Figure  13.  Addition 
of  carboxylic  acid  groups  to  bipyridine  ligand  of  dye  1  creating  dye  2,  results  in  the  significant  lowering  of 
the  orbital  energies  of  LUMO  -  LUMO+2  states.  Replacement  of  the  two  CN'  ligands  in  dye  2  by  an 
additional  bipyridine  ligand  creating  dye  3  results  in  a  set  of  three  of  doubly  degenerate  virtual  orbitals 
(LUMO  -  LUMO+5).  ' 

Results  of  the  IET  simulations  are  summarized  in  Figures  14-16.  In  order  for  the  IET  to  be 
competitive  with  the  ultrafast  intersystem  crossing  into  a  low-lying  5T  state  in  these  complexes,  the  IET 
characteristic  time,  x,  should  be  smaller  or  equal  to  100  fs.  ’  A  horizontal  red  line  in  Figures  14-16 
denotes  the  characteristic  time  x  =  100  fs.  Therefore,  only  the  initial  states  with  this  or  faster  rates  (i.e., 
lying  below  the  red  line)  will  be  capable  of  the  electron  injection  into  the  TiC>2  semiconductor. 

As  can  be  seen  in  Figure  14,  excited  dye  1  will  not  undergo  IET  efficiently,  as  the  majority  of  the 
characteristic  times  are  between  483.0  fs  -  17.7  ps.  The  only  exception  is  the  IET  from  LUMO+1  orbital 
in  the  bpy  perpendicular  attachment.  In  the  LUMO+1  case  the  characteristic  times  range  between  137.3  - 
101.6  fs,  suggesting  that  while  the  IET  is  not  very  efficient  in  1-Ti02  assemblies,  it  does  occur  for  certain 
excitations.  This  is  consistent  with  the  work  of  Meyer14  which  indicates  that  the  excited  state  charge 
transfer  in  1-Ti02  assemblies  occurs  via  a  combination  of  the  direct  sensitization  mechanism  (i.e.,  direct 
charge-transfer  excitation  between  the  dye  and  the  semiconductor)  as  well  as  the  indirect  mechanism  (i.e., 
IET).  It  is  also  interesting  to  note  that  the  IET  is  overall  more  efficient  when  the  plane  defined  by  the 
single  bipyridine  ligand  is  oriented  perpendicular  with  respect  to  the  Ti02  surface.  Additionally,  variation 
of  the  CN'  group  distance  from  the  TiCL  surface  between  2.12  -  2.33  A  does  not  significantly  impact  the 
calculated  IET  rates.  Therefore,  only  Ti-N  distance  of  2.33  A  was  considered  in  all  subsequent  IET 
simulations. 

Figure  15  summarizes  the  IET  rates  for  2-TiC>2  assemblies.  Both  bidentate  and  monodentate 
attachments  via  the  carboxylic  acid  anchoring  group  were  investigated,  as  well  as  the  attachment  via  the 
CN'  ligand.  There  are  several  conclusions  that  can  be  drawn  from  these  results.  First,  IET  via  the  CN' 
ligand  for  2-TiCL  assemblies  is  much  less  efficient  than  in  the  case  of  the  l-TiCL.  This  is  likely  due  to  the 
stabilization  of  the  MLCT  states  localized  on  the  bpy  ligand  caused  by  addition  of  the  carboxylic  acid 
groups  (see  Figure  13).  As  a  result,  IET  via  the  CN'  anchoring  group  will  not  be  competitive  with  the 
ultrafast  intersystem  crossing  events  in  these  complexes.  Second,  this  complex  displays  a  band-selective 
IET  when  attached  to  the  TiCL  via  the  carboxylic  acid  anchoring  group  in  both  monodentate  and  bidentate 
binding  modes.  Interestingly,  IET  via  the  monodentate  binding  mode  is  overall  more  efficient  than  via  the 
bindentate  binding  mode.  Finally,  comparison  of  the  IET  characteristic  times  among  l-TiCL  and  2-TiCL 
leads  us  to  conclude  that  indirect  sensitization  will  be  more  efficient  in  2-TiC>2  thanks  to  the  presence  of 
the  carboxylic  acid  anchoring  groups.  On  the  other  hand,  attachment  of  this  complex  via  the  carboxylic 
acid  group  will  eliminate  the  direct  sensitization  mechanism  that  occurs  when  the  complex  is  attached  to 
UO2  via  the  CN'  group.14 
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Characteristic  IET  times  for  3-TiC>2  assemblies  are  shown  in  Figure  16.  As  can  be  seen  in  Figure 
16,  dye  3  undergoes  a  band-selective  IET  for  all  attachments  investigated.  Interestingly,  attachment  via  the 
carboxylic  acid  anchoring  group  in  monodentate  binding  modes  results  in  the  fastest  IET  rates  on  average. 
Note  that  while  LUMO  -  LUMO+4  belong  to  the  Jt*  ligand- localized  states,  LUMO+15  (which  matches 
Kohn-Sham  orbital  LUMO+12)  corresponds  to  a  metal-centered  state  (see  Figure  8).  This  accounts  for  an 
ultrafast  IET  from  this  state  (x  =  12.4  fs)  for  CN'  attachment,  while  significantly  slower  IET  rates  (x  = 
302.9  fs  -  7.3  ps)  were  observed  for  the  3-TiC>2  assemblies  employing  carboxylic  acid  as  the  anchoring 
group. 

Interestingly,  band-selective  sensitization  was  observed  in  dye-TiCL  assemblies  for  all  three  dyes 
investigated.  While  our  results  do  not  elucidate  this  phenomenon  for  the  l-TiCL  system,  they  do  provide  an 
explanation  in  case  of  the  2-TiC>2  and  3-TiC>2  assemblies.  Band  1  (the  lowest  energy  band  in  the  UV-vis 
spectrum)  corresponds  to  the  excitations  into  the  LUMO  of  2  and  a  nearly  degenerate  set  of  LUMO  and 
LUMO+1  orbitals  of  3  (see  Figures  7  and  8).  These  orbitals  lie  at  the  edge  of  the  TiCL  conduction  band 
(see  Figure  13)  where  the  density  of  TiCL  states  is  low,  which  results  (1)  into  a  small  number  of  TiCL 
acceptor  states  available  to  couple  with  the  donor  states  (LUMO,  LUMO+1)  of  2  and  3  and  (2)  small 
driving  force  for  the  IET.  Our  results  are  consistent  with  the  experimental  observations  of  Ferrere  and 
coworkes  who  first  identified  this  band-selective  behavior  in  3-Ti02  assemblies.13  Moreover,  they  also 
provide  a  pathway  to  improving  the  efficiency  of  these  chromophores:  In  order  to  improve  the  IET 
efficiency,  one  needs  to  increase  the  energy  of  the  LUMO  and  LUMO+1  states  which  should  increase  the 
driving  force  for  the  electron  injection  as  well  as  provide  a  larger  number  of  the  TiCL  acceptor  states 
capable  of  coupling  with  the  dye  donor  states. 

6.  Conclusion 

In  this  work,  we  investigated  the  applicability  of  the  B3LYP  functional  to  accurate  prediction  of 
ground  state  multiplicity  of  pseudo-octahedral  iron  complexes.  Iron  complexes  are  of  great  interest  due  to 
potential  applications  as  molecular  switches,  data  storage  materials,  and  chromophores  in  dye-sensitized 
solar  cells.  Their  properties  intrinsically  depend  on  their  ground  state,  which,  depending  on  the  ligand 
character,  can  be  either  low-spin  or  high-spin.  In  general,  density  functional  theory  has  difficulties  with 
predicting  the  correct  ground  state  multiplicity  of  such  compounds,  as  pure  functionals  favor  the  low-spin 
states  and  hybrid  functionals  with  high  fraction  of  exact  Hartree-Fock  exchange  tend  to  favor  high-spin 
ground  states. 

We  confirmed  the  linear  relationship  between  the  high-spin/low-spin  energy  splitting,  A EHS/ls 
(AEhs/ls  =  Ehigh.spin-Elow_spin),  as  well  as  A Hhs/ls  and  A GHs/ls,  and  the  amount  of  exact  exchange  in  the 
B3LYP  functional  as  previously  observed  by  Reiher  and  coworkers.25  For  coordination  complexes  of  iron 
with  mostly  covalent  character,  such  as  Fe(II)-polypyridines,  A EHs/ls  displays  a  strong  dependence  on  the 
amount  of  exact  exchange  in  the  functional  and  B3LYP  tends  to  favor  high-spin  ground  states.  On  the 
other  hand,  A EHS/LS  for  ionic  complexes  displays  a  weaker  dependence  on  ch  and  B3LYP  artificially 
stabilizes  the  low-spin  states. 

We  confirmed  that  in  most  cases  investigated,  dependence  of  A EHS/ls  on  c/  is  linear  over  the  range 
of  0-25%  of  exact  exchange  admixture  in  B3LYP.  We  found  that  the  slope  of  this  dependence  strongly 
correlates  with  the  average  change  in  the  metal-ligand  bond  lengths  between  the  low-spin  and  high-spin 
states.  Moreover,  the  error  in  A EHs/ls  for  a  group  of  structurally  related  complexes  (i.e.,  complexes  that 
undergo  similar  change  in  the  metal-ligand  bond  lengths  between  the  high-spin  and  low-spin  states)  is 
systematic  and  the  calculated  values  of  A EHS/ls  correctly  reflect  the  trend  in  the  relative  stabilities  of  their 
low-spin  and  high-spin  states.  Inclusion  of  thermal  and  entropic  corrections  tends  to  stabilize  the  high-spin 
states,  but  nevertheless  results  in  the  identical  behavior.  This  systematic  behavior  allows  for  a  ground 
state  determination  of  an  arbitrary  pseudo-octahedral  complex  of  iron  by  comparing  the  calculated 
energy  differences  between  the  singlet  and  quintet  electronic  states  of  an  unknown  complex  to  the 
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energy  differences  of  structurally  related  complexes  with  a  known,  experimentally  determined 
ground  state. 

We  have,  furthermore,  studied  initial  excitation  and  subsequent  IET  among  three  Fe(II)- 
polypyridine  dyes  and  TiCF  anatase  semiconductor.  We  find  that  all  three  complexes  investigated  can 
undergo  IET  at  the  time  scale  competitive  with  the  ultrafast  intersystem  crossing  of  the  initially  excited 
'MLCT  states  into  the  low-lying  metal-centered  states  of  higher  multiplicities.  The  sensitization  is  more 
efficient  when  the  dyes  are  attached  to  surface  via  the  carboxylic  acid  anchoring  group  rather  than  CN". 
Attachment  via  the  carboxylic  acid  in  a  monodentate  binding  mode  also  results  in  more  efficient  IET  rates 
in  comparison  to  the  bidentate  attachment. 

All  three  complexes  investigated  display  band-selective  sensitization,  which  is  consistent  with 
previous  experimental  findings  by  Ferrere  and  coworkers.11  13  The  origin  of  the  band-selective 
behavior  is  attributed  to  the  poor  alignment  of  the  lowest  energy  excited  states  with  the  conduction 
band  of  the  Ti02  semiconductor.  This  results  in  a  poor  driving  force  for  the  injection  and  insufficient 
density  of  the  semiconductor  acceptor  states  available  for  coupling  with  the  dye  donor  states.  One  way  to 
improve  the  efficiency  of  Fe(II)-polypyridine  based  dyes  would  be  to  increase  the  energy  of  their  LUMO 
and  LUMO+1  states  with  respect  to  the  conduction  band  of  the  semiconductor.  This  could  be,  for  example, 
achieved  by  targeted  functionalization  of  their  polypyridine  scaffolds  by  electron  donating  groups. 

In  conclusion,  we  have  shown  that  DFT  is  a  valid  tool  that  can  be  employed  in  computational 
studies  of  Fe-based  coordination  compounds.  DFT,  TD-DFT,  and  quantum  dynamics  simulations  were 
also  successful  in  describing  and  elucidating  the  band-selective  IET  behavior  in  Fe(II)-polypyridine-Ti02 
nanoparticle  assemblies. 
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8.  Tables 


Table  1.  Experimental  and  calculated  spin  multiplicity  of  pseudo-octahedral  first  row  transition  metal 
complexes  investigated  in  this  work.  *  denotes  thermal  spin  crossover  complexes. 


Complex 

Ground  State 
(Experimental) 

Ground  State 
(A E,  B3LYP/SDD, 
6-311+G*) 

Ground  State  at 
273.15  K 

(AG,  B3LYP/SDD, 
6-311+G*) 

rFe(bpy)2(CN)2l° 

Singlet65 

Singlet 

Quintet 

[Fe(bpy)(CN)4l2' 

Singlet65 

Singlet 

Singlet 

[Fe(CN)6]4' 

Singlet66 

Singlet 

Singlet 

[Fe(bpy)2(NCS)2lu 

Singlet*67 

Quintet 

Quintet 

[Fe(bpy)2Cl2l<) 

Quintet68 

Quintet 

Quintet 

rFe(bpy)3f+ 

Singlet66 

Singlet 

Quintet 

[Fe(phen)3l2+ 

Singlet66 

Quintet 

Quintet 

rFe(Lt)2l2+ 

Singlet69 

Quintet 

Quintet 

[Fe(Me-Lt)212+ 

Quintet69 

Quintet 

Quintet 

[Fe(tpen)]2+ 

Singlet*70 

Quintet 

Quintet 

[Fe(terpy)2l2+ 

Singlet71 

Quintet 

Quintet 

[Fe(tren(py)3)l2+ 

Singlet72 

Quintet 

Quintet 

[Fe(tren(6-Me-py)3)]2+ 

Quintet72 

Quintet 

Quintet 

rFe(H20)6l2+ 

Quintet66 

Quintet 

Quintet 

[Fe(NH3)6l2+ 

Quintet73 

Quintet 

Quintet 

[Fe(NH3)6l3+ 

Sextet73 

Doublet 

Sextet 

rFeF6l3' 

Sextet74 

Sextet 

Sextet 

rFeClfil3' 

Sextet75 

Sextet 

Sextet 

[Fe(pyrol)3tren]° 

Doublet76 

Doublet 

Sextet 

[Mn(pyrol)3tren]° 

Quintet77 

Quintet 

Quintet 

rCo(bpy)3l2+ 

Quartet78 

Quartet 

Quartet 

rCo(terpy)2l2+ 

Doublet*79 

Quartet 

Quartet 

rcoF6f- 

Quintet80 

Quintet 

Quintet 

rCo(CN)6l3- 

Singlet81 

Singlet 

Singlet 

rCo(NH3)6l3+ 

Singlet82 

Singlet 

Singlet 
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Table  2.  Table  of  average  change  in  metal  -  ligand  bond  lengths  AR( metal- ligand)  in  Angstroms  for  low 
spin  to  high  spin  transition  versus  slope  of  the  scan  over  the  exact  exchange  for  various  pseudo-octahedral 
complexes.  *  [0s(bpy)3]2+  quintet  changes  electronic  states  over  the  exact  exchange  admixture 
investigated. 


AR 

Slope  of  Scan 

Slope  of  Scan 

Slope  of  Scan 

Complex 

(metal- 

d(AE) 

cl  (AH) 

d(AG) 

ligand)  [A] 

dcx 

dcx 

dcx 

[Fe(CN)6]4' 

0.3336 

-214.17 

-212.51 

-208.92 

[Fe(bpy)(CN)4]2' 

0.2726 

-185.56 

-184.48 

-181.69 

[Fe(bpy)2(CN)2l° 

0.2337 

-161.42 

-160.73 

160.02 

[Co(CN)6]3- 

0.2295 

-145.62 

-143.93 

-137.59 

[Fe(Lt)2r 

0.2266 

-160.18 

-158.10 

-161.41 

[Fe(tren(py)3)]2+ 

0.2265 

-147.79 

-146.01 

-141.98 

[Fe(terpy)2]2+ 

0.2149 

-156.18 

-154.86 

-152.45 

[Fe(tren(6-Me-py)3)]2+ 

0.2125 

-142.97 

-141.98 

-142.88 

[Fe(tpen)]2+ 

0.2044 

-125.25 

-124.52 

-117.40 

[Fe(bpy)3]2+ 

0.2013 

-143.56 

-143.20 

-141.62 

[Fe(bpy)2(NCS)2]w 

0.2003 

-125.76 

-124.30 

-124.48 

[Fe(phen)3]2+ 

0.1994 

-139.34 

-139.67 

-138.29 

[Fe(Me-Lt)2|2+ 

0.1991 

-142.16 

-149.73 

-156.79 

[Fe(NH3)6l2+ 

0.1827 

-91.95 

-92.05 

-90.70 

[Fe(bpy)2Cl2|u 

0.1755 

-126.90 

-127.05 

-127.54 

[Fe(pyrol)3tren]° 

0.1655 

-97.80 

-96.84 

-92.48 

[Co(NH3)6]3+ 

0.1561 

-88.39 

-90.37 

-93.56 

[Fe(NH3)6l3+ 

0.1498 

-89.33 

-92.26 

-93.07 

[Co(bpy)3l2+ 

0.1290 

-89.33 

-89.55 

-87.59 

[Fe(H20)6]2+ 

0.1230 

-62.12 

-62.24 

-63.08 

[Co(terpy)2]2+ 

0.1136 

-93.29 

-90.49 

-97.06 

[Mn(pyrol)3tren]° 

0.1095 

-58.44 

-58.88 

-57.12 

[FeCl6l3- 

0.0908 

-54.97 

-54.51 

-49.97 

[CoF6]3' 

0.0638 

-45.42 

-38.17 

-30.97 

[FeF6]3' 

0.0610 

-48.25 

-47.47 

-42.93 

[Ru(bpy)3]2+ 

0.2760 

-109.44 

-108.74 

-112.39 

[Os(bpy)3]2+ 

0.1372 

-16.74* 

-9.21* 

-17.89* 
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Table  3.  Table  of  metal  ligand  bond  lengths  for  the  complexes  investigated,  coordination  site  labeling 
corresponds  with  Figure  6. 


[Fe(bpy)(CN)4]2- 

[Fe(bpy-dca)(CN)4]2' 

[Fe(bpy-dca)2(CN)2] 

Iron-Ligand  Bond 

1 

2 

3 

Length  (A) 

Length  (A) 

Length  (A) 

Fe-Ll  (CN  ) 

1.959 

1.955 

1.942 

Fe-L2  (CN') 

1.988 

1.977 

1.942 

Fe-L3  (CN',  bpy) 

1.959 

1.955 

2.020 

Fe-L4  (CN',  bpy) 

1.989 

1.977 

1.996 

Fe-L5  (bpy) 

1.992 

1.959 

2.020 

Fe-L6  (bpy) 

1.992 

1.963 

1.998 
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9.  Figure  Captions 

Figure  1.  Fe(II)  complexes  investigated  in  this  work 
Figure  2.  Heteropyridine  ligands  considered  in  this  study. 


Figure  3.  Dependence  of  the  energy  difference,  AEHSILS  =  Ehj  h  in  -  Elow  in  (kcal/mol),  on  the  exact 

exchange  admixture  (c,)  in  the  B3LYP  functional  for  iron  pseudo-octahedral  complexes  with  various 
ligands:  [Fe(CN)6]4-  (1),  [Fe(bpy)2(CN)2]°  (2),  [Fe(bpy)3]2+  (3),  [Fe(NH3)6]3+  (4),  [FeCl6]3-  (5).  A EHS/LS>  0 
corresponds  to  the  low-spin  ground  state  (singlet  or  doublet),  A EHS/ls<  0  corresponds  to  the  high-spin 
ground  state  (quintet  or  sextet). 

Figure  4.  Plot  of  average  change  in  metal  to  ligand  bond  lengths  in  transition  from  low-spin  to  high-spin 
states  versus  slope  of  the  scan  over  exact  exchange  for  data  shown  in  Table  2  along  with  a  plot  of  the 
linear  regression,  R2=0.93. 


Figure  5.  Dependence  of  the  energy  difference,  AEHS/LS  =  Ehigh_spin  -  Elow_spin  AEHS/LS  =  Ehigh_spin  - 

Eiow -spin  (kcal/mol),  on  the  exact  exchange  admixture  (c/)  in  the  B3LYP  functional  for  several  pseudo- 
octahedral  iron(II)-polypyridine  complexes:  [Fe(bpy)2(CN)2]°  (1),  [Fe(bpy)3]2+  (2),  [Fe(tren(py)3)]2+  (3), 
[FefbpyHNCSy0  (4),  [Fe(tren(6-Me-py)3)]2+  (5),  and  [Fe(bpy)2Cl2]0  (6).  AEHs/ls  >  0  A EHS/LS  >  0 
corresponds  to  the  singlet  ground  state,  A EHs/ls  <  0  A EHS/LS  <  0  corresponds  to  the  quintet  ground  state.  A 
linear  regression  is  plotted  for  each  complex  based  on  the  constant  slope  formula  A EHS/LS  =  I  - 140.2^  . 

Figure  6.  Optimized  structures  of  dyes  1,2,  and  3  with  each  coordination  site  labeled  (L1-L6). 

Figure  7.  Simulated  absorption  spectra  for  (from  top  to  bottom)  dyes  1,  2,  and  3  with  Lorentzian 
broadening  of  HWMH  =  0.12  eV  with  important  excitations  (fosc  >  0.01)  in  the  visible  region  labeled, 
calculated  with  B3LYP  TD-DFT  in  PCM  (acetonitrile). 

Figure  8.  Relevant  particle  states  for  the  major  excitations  (fosc  >  0.01)  from  the  calculated  spectra  shown 
in  Figure  7.  Kohn-Sham  orbitals  are  classified  by  the  absorption  band  (columns)  and  by  similarity  in  nodal 
structure  (rows). 

Figure  9.  Optimized  structures  of  anatase  (101)  with  pyridine-4-carboxylic  acid  binding  via  a 
monodentate  (a)  and  bidentate  (b)  binding  mode,  and  hydrogen  cyanide  (c). 

Figure  10.  Final  slab  models  of  dye-nanoparticle  assemblies  consisting  of  [Fe(bpy-dca)2(CN)2]  attached  to 
anatase  (101)  via  a  monodentate  carboxylic  acid  binding  mode  in  the  following  nonequivalent 
orientations:  bpy  parallel  (left)  and  bpy  perpendicular  (right). 

Figure  11.  Structure  of  1-Ti02  assembly  employing  CN'  as  the  anchoring  group. 

Figure  12.  Density  of  states  for  1-Ti02  anatase  (101)  slab  model  (left)  and  enlarged  conduction  band 
(right)  show  the  following:  total  density  of  states  (blue  line),  projected  density  of  states  on  the  dye  (black 
line),  and  the  energy  levels  of  the  dye  in  vacuum  using  extended  Hiickel  (black  level  set  lines).  Gaussian 
line-shape  (HWHM  =  0.05  eV)  used  for  convolution. 
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Figure  13.  Discrete  energy  levels  of  the  dyes  (left)  in  vacuum,  using  extended  Hiickel  level  of  theory, 
showing  the  evolution  of  the  similar  molecular  orbitals,  going  from  the  smallest  to  the  largest  dye.  The 
DOS/pDOS  for  the  [Fe(bpy-dca)2(CN)2]/Ti02  nanoparticle  system  is  shown  for  reference  (right). 

Figure  14.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of  [Fe(bpy)(CN)4]2' 
attached  to  anatase  (101)  slab  in  both  the  “bpy  parallel”  and  “bpy  perpendicular”  orientations  via 
monodentate  cyanide  anchoring  groups  with  varying  nitrogen  titanium  distances  (2.33,  2.22,  and  2.12  A), 
determined  by  either  exponential  or  biexponential  fitting. 

Figure  15.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of  [Fe(bpy- 
dca)(CN)4]2"  attached  to  anatase  (101)  slab  by  monodentate  carboxylic  acid,  bidentate  carboxylic  acid,  and 
cyanide  anchoring  groups  (in  both  the  “bpy  parallel”  and  “bpy  perpendicular”  orientations),  determined  by 
either  exponential  or  biexponential  fitting. 

Figure  16.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of 
[Fe(bpy-dca)2(CN)2]  attached  to  anatase  (101)  slab  by  monodentate  carboxylic  acid  (in  both  the  “bpy 
parallel”  and  “bpy  perpendicular”  orientations),  bidentate  carboxylic  acid  (in  both  the  “bpy  parallel”  and 
“bpy  perpendicular”  orientations),  and  cyanide  anchoring  groups,  determined  by  either  exponential  or 
biexponential  fitting. 
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10.  Figures 


Figure  1.  Fe(II)  complexes  investigated  in  this  work 
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Figure  3.  Dependence  of  the  energy  difference,  A EHSILS  =  Ehjgh_spin- Elow_spin  (kcal/mol),  on  the  exact 

exchange  admixture  (c,)  in  the  B3LYP  functional  for  iron  pseudo-octahedral  complexes  with  various 
ligands:  [Fe(CN)6]4-  (1),  [Fe(bpy)2(CN)2]°  (2),  [Fe(bpy)3]2+  (3),  [Fe(NH3)6]3+  (4),  [FeCl6f  (5).  A EHS/LS>  0 
corresponds  to  the  low-spin  ground  state  (singlet  or  doublet),  A EHS/LS<  0  corresponds  to  the  high-spin 
ground  state  (quintet  or  sextet). 
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0.35 


♦  Fe,  Co,  Mn  complexes  ♦ 


Slope  of  the  Scan 

Figure  4.  Plot  of  average  change  in  metal  to  ligand  bond  lengths  in  transition  from  low-spin  to  high-spin 
states  versus  slope  of  the  scan  over  exact  exchange  for  data  shown  in  Table  2  along  with  a  plot  of  the 
linear  regression,  R2=0.93. 
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pure  DFT  limit  HF  limit - ►  cn=1 

Figure  5.  Dependence  of  the  energy  difference,  AEHSI!S  =  Ehigh_spin  -  Elow_spin  A EHS/LS  =  Ehigh_spin  - 

E ioW -spin  (kcal/mol),  on  the  exact  exchange  admixture  (cj)  in  the  B3LYP  functional  for  several  pseudo- 
octahedral  iron(II)-polypyridine  complexes:  [Fe(bpy)2(CN)2]°  (1),  [Fe(bpy)3]2+  (2),  [Fe(tren(py)3)]2+  (3), 
[Fe(bpy)2(NCS)2]0  (4),  [Fe(tren(6-Me-py)3)]2+  (5),  and  [Fe(bpy)2Cl2]°  (6).  A EHS/ls  >  0  A EHS/LS  >  0 
corresponds  to  the  singlet  ground  state,  AEHS/ls  <  0  A EHS/LS  <  0  corresponds  to  the  quintet  ground  state.  A 

linear  regression  is  plotted  for  each  complex  based  on  the  constant  slope  formula  A EHS/LS  = /-140.2c, 
AEhs/ls  =  I  ~  140.2c! 
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Figure  7.  Simulated  absorption  spectra  for  (from  top  to  bottom)  dyes  1,  2,  and  3  with  Lorentzian 
broadening  of  HWMH  =  0.12  eV  with  important  excitations  (fosc  >  0.01)  in  the  visible  region  labeled, 
calculated  with  B3LYP  TD-DFT  in  PCM  (acetonitrile). 
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Figure  8.  Relevant  particle  states  for  the  major  excitations  (fosc  >  0.01)  from  the  calculated  spectra  shown 
in  Figure  7.  Kohn-Sham  orbitals  are  classified  by  the  absorption  band  (columns)  and  by  similarity  in  nodal 
structure  (rows). 
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Figure  9.  Optimized  structures  of  anatase  (101)  with  pyridine-4-carboxylic  acid  binding  via  a 
monodentate  (a)  and  bidentate  (b)  binding  mode,  and  hydrogen  cyanide  (c). 
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Figure  10.  Final  slab  models  of  dye-nanoparticle  assemblies  consisting  of  [Fe(bpy-dca)2(CN)2]  attached  to 
anatase  (101)  via  a  monodentate  carboxylic  acid  binding  mode  in  the  following  nonequivalent 
orientations:  bpy  parallel  (left)  and  bpy  perpendicular  (right). 
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N-Ti:  used  2.3,  2.2,  and  2.1  A 


Figure  11.  Structure  of  1-Ti02  assembly  employing  CN'  as  the  anchoring  group. 
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Figure  12.  Density  of  states  for  1-Ti02  anatase  (101)  slab  model  (left)  and  enlarged  conduction  band 
(right)  show  the  following:  total  density  of  states  (blue  line),  projected  density  of  states  on  the  dye  (black 
line),  and  the  energy  levels  of  the  dye  in  vacuum  using  extended  Hiickel  (black  level  set  lines).  Gaussian 
line-shape  (HWHM  =  0.05  eV)  used  for  convolution. 
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Figure  13.  Discrete  energy  levels  of  the  dyes  (left)  in  vacuum,  using  extended  Hiickel  level  of  theory, 
showing  the  evolution  of  the  similar  molecular  orbitals,  going  from  the  smallest  to  the  largest  dye.  The 
DOS/pDOS  for  the  [Fe(bpy-dca)2(CN)2]/Ti02  nanoparticle  system  is  shown  for  reference  (right). 
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Figure  14.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of  [Fe(bpy)(CN)4]2' 
attached  to  anatase  (101)  slab  in  both  the  “bpy  parallel”  and  “bpy  perpendicular”  orientations  via 
monodentate  cyanide  anchoring  groups  with  varying  nitrogen  titanium  distances  (2.33,  2.22,  and  2.12  A), 
determined  by  either  exponential  or  biexponential  fitting. 
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Figure  15.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of  [Fe(bpy- 
dca)(CN)4]2"  attached  to  anatase  (101)  slab  by  monodentate  carboxylic  acid,  bidentate  carboxylic  acid,  and 
cyanide  anchoring  groups  (in  both  the  “bpy  parallel”  and  “bpy  perpendicular”  orientations),  determined  by 
either  exponential  or  biexponential  fitting. 


38 


UJ 

o 

’-I— < 

<S) 

CD 

-i— < 

o 

cc 

cC 

JZ 

o 


105 

104 

103 

102 

101 

10° 


Monodentate  CA,  Parallel 
Monodentate  CA,  Perpendicular 
Bidentate  CA,  Parallel 
Bidentate  CA,  Perpendicular 

CN 


. 


Figure  16.  Characteristic  interfacial  electron  transfer  times  for  relevant  particle  states  of 
[Fe(bpy-dca)2(CN)2]  attached  to  anatase  (101)  slab  by  monodentate  carboxylic  acid  (in  both  the  “bpy 
parallel”  and  “bpy  perpendicular”  orientations),  bidentate  carboxylic  acid  (in  both  the  “bpy  parallel”  and 
“bpy  perpendicular”  orientations),  and  cyanide  anchoring  groups,  determined  by  either  exponential  or 
biexponential  fitting. 
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